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I.  INTRODUCTION 


In  recent  years,  the  Chemical  Kinetics  Department  of  The  Aerospace 

Corporation  Aerophysics  Laboratory  has  prepared  a series  of  reports  that 

1-3 

review  the  kinetics  of  hydrogen  halide  lasers  HF,  DF,  and  HCl.  The 
most  recent  of  these  devoted  to  the  HF  system  was  published  in  April  1976 
and  covered  the  literature  through  June  1975.  Most  of  the  1976  report  and 
its  recommendations  still  appear  to  be  valid,  although  kinetic  work  continues 
to  be  done  on  the  HF  system,  particularly  in  regard  to  energy -transfer 
processes.  Therefore,  rather  than  attempting  a complete  revision  of  that 
1976  report,  we  present  a supplement  that  includes  discussion  of  work  that 
has  appeared  during  the  past  two  years,  with  revised  conclusions  where  ap- 
propriate. A list  of  recommended  rate  coefficients,  based  on  all  the  litera- 
ture through  September  1977,  appears  in  the  Appendix. 


II.  DISSOCIATION -RECOMBINATION  REACTIONS 


Recent  measurements  by  Ultee  of  fluorine  atom  recombination  confirm 
previous  indications  that  the  rate  coefficient  is  srriall  compared  with  those 
for  other  atomic  species.  Ultee  used  electron  spin  resonance  (esr)  to  mea- 
sure fluorine  atom  decay  at  room  temperature  down  a quartz  fast-flow  system 
in  which  F atoms  were  produced  by  microwave  discharge  through  a F^- 
carrier  gas  mixture.  Data  were  reported  at  total  pressures  of  8 to  34  Torr 
in  mixtures  of  F_  and  He  that  consisted  of  2 to  10%  F,.  A value  of 
k = 2.2±0.3X  10  /mole^- sec  was  reported;  preliminary  data  indicated 

< k ^ < k This  is  reasonably  consistent  with  the  value  of  k^^  = 

1 X 10^"^  cm  /mole  -sec  recommended  by  Cohen  and  Bott.  ^ However,  in 

He  .A.  r 

Ref.  Z it  was  assumed  k = k . In  view  of  Ultee's  results,  perhaps  it 

H e r • 

would  be  more  correct  to  assume  that  k = 2k  * 

Previous  reports  in  this  series  have  not  dealt  with  wall  recombination 
of  atomic  species.  In  recent  years,  several  authors  have  reported  measure- 
ments of  F-atom  combination  on  surfaces.  A summary  of  the  existing  litera- 
ture on  this  subject  is  given  in  Table  1.  In  general,  results  indicate  a re- 

-5  -4 

combination  coefficient  of  10  - 10  for  pyrex,  coated  or  uncoated  quartz, 

and  alumina.  Metallic  surfaces  are  somewhat  more  efficient,  with  ys  rang- 
ing from  10'^  to  10“^.  In  several  cases,  measurements  were  made  over 
a sufficient  temperature  range  to  report  Arrhenius  parameters  for  the  re- 
combination rate  coefficient.  For  quartz,  molybdenum  glass,  and  gold,  there 
is  a slight  positive  activation  energy  of  about  3 kcal/mole;  for  copper,  there 
is  a negative  activation  energy  of  about  5 kcal/mole.  Heterogeneous  atom 
disappearance  can  be  the  result  of  several  processes,  including  simple  diatom 
formation  catalyzed  by  the  surface  and  chemisorption  of  atoms  on  the  surface, 
which  results  in  fluorides  of  the  wall  material.  It  is  impossible  to  deduce  from 
the  Arrhenius  parameters  which  type  of  process  is  occurring.  Either  simple 
catalysis  or  wall  chemisorption  can  have  a positive  or  negative  activation 
energy. 
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Table  1.  Heterogeneous  F Atom  Loss  Probabilitie 
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There  have  been  no  important  recent  results  in  or  HF  dissociation 

recombination.  Surface  catalyzed  recombination  of  H atoms  was  reviewed 

12 

in  1973  by  Jones,  MacKnight,  and  Teng. 
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III.  PUMPING  REACTIONS 


Interest  has  grown  recently  in  the  kinetics  of  fluorine  reactions,  and  an 
increasing  number  of  experimental  studies  have  been  devoted  to  such  pro- 
cesses. The  difficulties  of  measuring  the  rate  of  the  F + H2  reaction  were 

discussed  in  Ref.  1.  Since  that  review  was  compiled,  Goldberg  and 
13 

Schneider  have  reported  a measurement  of  for  the  process 

F + - HF  + H (4)* 

They  studied  the  reaction  using  a fast-flow  reactor  at  298  K with  esr  detec- 
tion to  measure  F-atom  concentration  as  a function  of  distance  along  the 
2.2  cm  i.d.  quartz  flow  tube.  A microwave  discharge  through  F2  was  the 
source  of  F atoms.  Calibration  indicated  that  the  F2  was  nearly  100%  dis- 
sociated, which  avoided  any  of  the  complications  that  would  result  from  F^ 

13  3 ^ 

being  present.  A value  of  k.  = 1.0  • 0.08  x 10  cm  /mole-sec  was  obtained. 

^ 13  3 

This  value  compares  favorably  with  the  value  of  1.5  x 10  cm  /mole-sec 
recommended  by  Cohen  and  Bott.  In  the  same  facility,  but  under  conditions 
in  which  the  F2  was  not  completely  dissociated,  Goldberg  and  Schneider 
measured  the  value  of  k^  for  the  process 

H + F2  - HF  + F (5) 

In  this  case,  the  system  was  somewhat  more  complicated,  resulting  in 

greater  uncertainty  in  the  results.  They  reported  a value  of  k^  = 2.6  • 0.6 

12  3 ^ 

X 10  cm  /mole-sec.  The  range  of  uncertainty  encompasses  the  value  of 

2.  1 X 10  reported  by  Albright  et  al.  and  recommended  by  Cohen  and  Bott. 


Reactions  here  are  numbered  to  be  consistent  with  Ref.  1. 


Manning  et  al.  studied  the  hydrogen  abstraction  reactions  by  atomic 

fluorine  from  a series  of  hydrocarbons  and  fluorinated  hydrocarbons  at 

1 8 

303  K.  Neutron  bombardment  was  used  to  produce  F from  the  stable  isotope 
19 

F.  Pairs  of  reagents  in  the  presence  of  an  excess  of  moderator  bath  were 

18 

reacted  with  F: 


4 0 4 0 

F + RH  - H ^F  + R 

18_  TP  r'  TP 

F + C3F^  - C3F^  F 

1 8 

The  quantity  of  H F,  measured  by  radiochemical  techniques,  was  the  only 
variable  monitored.  Analysis  of  the  data  required  knowledge  of  the  fraction 
of  hot  F atoms  that  were  cooled  to  thermal  temperatures  before  reaction; 
this  number,  a constant  for  a given  series  of  experiments,  was  found  to  be 
0.85  J^O.02.  A value  of  0.39  was  reported  for  the  ratio  of  rate  coefficients 
kp-j.fj2/^F+CH4*  Since  the  best  value  for  appears  to  be  4 1 

X lO^^cm^/mole-sec  (see  Ref.  1 for  discussion),  this  yields  a value  of 
k^  = 1.6  X lO^"^,  in  good  agreement  with  the  previously  recommended  value. 

The  room-temperature  relative  distribution  of  HF(v)  formed  in  reac- 
tion (4)  has  been  reasonably  well  established  after  several  years  of  con- 
tinual improvements  in  experimental  techniques.  Perry  and  Polanyi,  con- 
tinuing the  careful  work  carried  out  in  Polanyi *s  laboratory,  used  the 
arrested  relaxation  infrared  chemiluminescence  technique  to  obtain  relative 

rates  of  production  of  v = 0 . . . . 3 of  0 . 04:0.  28:1 . 0:0.  55,  in  good  agreement 

17 

with  the  recommendations  of  Ref.  1.  Berry  reported  a somewhat  dif- 
ferent distribution  of  0.06:0.29:1.0:0.63  by  using  the  quite  different  tech- 
nique of  measurements  of  the  relative  gain  coefficients  of  individual  HF 
laser  transitions.  He  suggested  that  the  discrepancy  in  results  for 

ov-k.,  might  in  some  cases  be  due  to  vibrational  relaxation  of 
4(v=3)  4(v=2)  ^ 
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1 8 

product  HF  (cf.  Coombe  and  Pimentel  ) or  differences  in  translational 
energies  of  the  reagents.  However,  Krogh,  Stone,  and  Pimental  ^ ex- 
amined in  detail  the  effects  of  vibrational  relaxation  on  the  results  of  Coombe 
and  Pimentel  and  concluded  that  energy-transfer  processes  did  not  interfere 
with  the  measurements.  Instead,  they  suggested  rotational  disequilibrium 
as  a possible  source  for  the  discrepancy. 


21 

Moehlmann  and  McDonald  obtained  very  similar  results  with  the 

arrested  relaxation  infrared  chemiluminescence  technique  similar  to  that  of 

Polanyi  and  co-workers.  They  reported  relative  population  rates  for 

V = 1 , 2,  3 to  be  0.  28:1. 0:0. 47.  In  view  of  the  near  agreement  among  all  the 

experiments  save  one,  it  seems  best  to  retain  the  recommendation  of  Ref.  1; 

that  t A /o\ / = 0-  17:0. 55:0.  28  = 0.31:1.0:0.51  at  room  temperature. 

4(1)  4(2)  4(3)  ^ 

No  measurements  of  the  temperature  dependences  of  these  ratios  have  been 
reported  since  the  work  of  Coombe  and  Pimentel. 


The  F + reaction  (and  its  isotopic  analogs)  has  been  the  subject  of 

^ 22 
several  theoretical  studies  in  recent  years.  Connor,  Jakubetz,  and  Manz 

analyzed  the  population  distributions  predicted  by  one-dimensional  trajectory 

calculations  using  information  theoretic  surprisal  plots.  Kormornicki, 

23 

George,  and  Morokuma  studied  the  effect  of  multiple  electronic  surfaces  on 

24 

the  reactivity  of  F + Rebentrost  and  Lester  examined  nonadiabatic 

effects  of  interaction  potentials  for  the  various  states  of  the  triatomic  inter- 

25 

mediate  in  the  reaction.  Schatz,  Bowman,  and  Kuppermann  presented  and 
compared  exact  quantum,  quasi -classical,  and  semiclassical  reaction  proba- 
bilities and  rate  coefficients  for  the  collinear  reaction.  Polanyi  and 
26 

Schreiber,  continuing  earlier  studies  by  Polanyi  and  co-workers,  reported 
an  extensive  series  of  investigations  of  the  detailed  dynamics  of  the  F + H2 
reaction  on  five  different  potential  surfaces;  they  examined  the  effect  of  re- 
reagent vibrational  or  translational  energy  on  the  product  energy  distribu- 
tions. All  of  these  studies  can  provide  useful  insights  into  the  fundamental 
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chemical  dynamics  of  a system  such  as  F + However,  inasmuch  as  they 

do  not  provide  any  information  over  and  above  the  experimental  data  neces- 
sary  for  evaluating  the  rate  parameters  needed  for  modeling  laser  systems, 
they  are  not  discussed  in  detail  here. 

27 

Polanyi,  Sloan,  and  Wanner  studied  the  H + reaction  using  the 

ar rested-relaxation  infrared  chemiluminescence  technique.  They  examined 

the  effect  on  product  energy  distribution  of  separately  enhanced  translational 

energy  (in  the  H atom)  and  vibrational-rotational  energy  (in  the  F^  molecule). 

When  the  translational  energy  of  the  H was  increased  (by  formation  of  the 

species  in  an  oven  at  temperatures  up  to  2800  K),  the  excess  energy  was 

channeled  fully  into  translational  and  rotational  energy  of  the  products;  the 

vibrational  distribution  was  altered  slightly,  with  the  peak  of  production 

vibrational  energy  shifted  from  HF(v  = 6)  to  HF(v  = 5).  Heating  the  F^  to 

900  K produced  no  measurable  effect  on  the  product  energy  distribution. 

Although  these  results  are  not  directly  comparable  to  experiments  in  which 

the  reagents  are  all  at  the  same  temperature,  they  do  suggest  that,  if  the 

reagents  are  thermally  equilibrated,  there  will  be  a decrease  in  energy 

deposited  in  HF  vibration  as  the  reagent  temperature  is  increased.  Calcu- 

2 8 

lations  by  Wilkins  predict  that  the  fraction  of  reagent  energy  channeled 
into  vibration  of  HF  will  be  almost  constant  from  300  to  1000  K. 
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IV.  ENERGY-TRANSFER  PROCESSES 


There  are  two  persistent  problem  areas  in  HF  vibrational  relaxation: 

(1)  relaxation  by  H atoms,  and  (Z)  self- relaxation  by  a combination  of  V-V 
and  V-R,  T energy  transfer.  These  processes  continue  to  demand  attention 
because  (1)  they  are  the  fastest  deactivation  processes  in  the  HF  system, 
and  therefore  are  principal  factors  in  the  loss  of  laser  power;  and  (2)  they 
are  the  most  difficult  to  resolve  experimentally,  particularly  when  upper 
vibrational  levels  of  HF  are  involved. 

A.  HF  SELF  RELAXATION 

29 

^ Recently,  Wilkins  completed  the  first  three-dimensional  trajectory 

study  of  HF-HF  interactions.  Because  his  studies  provide  a convenient 
unified  framework  in  which  to  consider  all  the  possible  HF-HF  processes, 
a summary  of  his  conclusions  is  presented  first.  The  experimental  data 
are  reviewed  and  discussed  here  in  the  context  of  the  theoretical  studies. 

An  important  conclusion  of  the  trajectory  studies  is  that  a favored 
process  in  the  collision  between  two  HF  molecules,  one  of  which  is  vibra- 
tionally  excited,  is  the  internal  conversion  of  vibrational-to-rotational  energy, 
with  negligible  energy  being  transferred  to  the  collision  partner.  For  a 
given  pair  of  molecules,  HF(Vj,  J^)  and  HF(v2>  J2)>  collision  products  are 
HF(vj-Avj,  Jp  and  HF(v2,  with  Av^  varying  from  1 to  v^  and  varying 

over  a half  dozen  states  from  the  highest  one  energetically  possible  on  down, 
all  with  approximately  equal  probability.  The  resulting  states  are  sum- 
marized in  Table  2 for  = 2,  the  most  probable  state  at  300  K.  The  un- 
certainties in  the  rate  coefficients  reflect  the  statistically  small  number  of 
trajectories  giving  the  indicated  product  states.  Within  the  uncertainties,  it 

appears  that  the  rate  coefficients  are  all  approximately  equal:  between  1 and 
12  3 

7 xlO  cm  /mole-sec  for  any  given  process.  The  rate  coefficients  are 
listed  in  Table  3 for  the  J-specific  coefficients  summed  over  all  J states  and 
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Table  2,  Theoretical  Vibrational-to- Rotational  Rate  Coefficients  and  Product 
Channels  for  HF(vi,  Jj)  + HF(v2,  J2)^HF(vi,  Ji)  + HF(v^.  J^)  at 
T = 300,  V2  = = 0,  Jj  = 2^ 
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after  taking  into  account  the  rotational  state  populations  at  equilibrium. 


Table  3.  Theoretical  Rate  Coefficients  kv4,v2;  vf,v^  for  Vibrational-to-Rotational 
Energy  Transfer  in  HF{vi)  + HF(v2)  Collisions  at  T = 300 
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weighted  according  to  the  Boltzmann  populations  of  rotational  levels  at  300  K. 
The  table  shows  that  a collision  between  HF(v)  and  HF(0)  will  give  HF(v^) 

+ HF(0),  where  v^  varies  from  v-1  to  0.  Single -quantum  transitions,  i.e., 


v^  = v-1,  are  the  most  probable.  Multiquantum  transitions  are  about  one- 
half  as  likely  as  single-quantum  transitions,  but  there  is  little  difference 


between  the  probabilities  of  losing  2,3.  . .v  quanta.  Similar  results  are 
obtained  when  both  HF  molecules  are  vibrationally  excited,  except  that  now 
there  is  a considerably  larger  number  of  possible  exit  channels.  One  or 
both  molecules  can  lose  vibrational  energy,  which  is  converted  intramolecu- 
larly  to  rotational  energy.  For  example,  the  vibrational  energy  lost  by  each 
molecule  is  largely  converted  to  rotational  energy  of  that  same  molecule. 

The  probability  of  both  molecules  losing  vibrational  energy  appears  to  be 
smaller  than  that  for  only  one  molecule  to  lose  vibrational  energy.  However, 
there  is  no  simple  correlation  between  the  change  in  v and  the  rate 
coefficients. 

Unlike  V-R,  T energy  transfer,  in  which  large  Av  processes  were 
observed  with  high  probability,  in  V-V  exchange,  the  trajectory  calculations 
predicted  that  only  Av  = 1 is  likely.  Furthermore,  the  probabilities  for 
HF(v)  + HF(v  0 HF(v  - 1)  + HF(v'  + 1)  were  found  to  decrease  with  increasing  v. 
Within  approximately  25%,  the  results  of  Wilkins  for  V-V  exchange  (up  to 

5jC 

V = 5)  can  be  summarized  by  the  expressions 


k 


1,  v;0, v+1 


^3.6x10 


(A) 


k 


V,  v;v-  1 , v+ 1 


“3X10 


v>  1 


(B) 


R.  L.  Wilkins,  private  communication. 
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In  processes  of  this  sort,  the  rotational  states  of  the  HF  molecules  change 
during  collision;  thus,  AE  for  any  process  is  considerably  smaller  than  what 
would  be  calculated  by  assuming  no  change  in  rotational  energy  across  the  re- 
action. (It  should  be  stressed  that  the  above  two  analytic  expressions  for  k 
represent  simple,  approximate  fits  to  the  results  of  the  trajectory  calcula- 
tions; the  V dependences  have  no  theoretical  significance.) 

29 

The  results  of  Wilkins  for  V-V  and  V-R,  T transfer  ' taken  together 
suggest  that,  for  small  v,  both  kinds  of  processes  are  approximately  equally 
probable,  contrary  to  earlier  assumptions  that  V-V  was  considerably  faster 
than  V-T.  For  example,  in  collisions  between  HF(v=l)  and  HF(v=l),  Wilkins 
reported  the  rate  coefficients  given  in  Table  4 for  the  various  possible 
exit  channels. 


Table  4.  Calculated  Rate  Coefficients  for 
HF(  1 ) + HF(  1 ) - HF(vp  + HF(v2  ) 
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"^2 

. / / 
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0 

4.  1 ± 1.3  ) 

[ V-R,  T 

0 

0 

1.  5 ± 0.9  ) 

0 

2 

1. 2 ± 0.  2 V-V 

Furthermore,  whereas  the  overall  V-R,T  rates  increase  with  v,  the  V-V 
rates  decrease. 

In  experiments  designed  to  measure  the  V-V  rate  of  HF  relaxation 
that  have  been  conducted  to  date,  one  observes  the  decay  of  (HF(v)],  which 


is  due  to  two  competing  processes,  a V-V  exchange  (endothermic)  and  a 
V-R,T  transfer  (exothermic) 


HF(v)  + HF(0)  - HF(v.Av).+  HF(1)  (V-V) 

HF(v,  J)  + HF(0)  - HF(v  J 0 + HF(0)  ( V-R,  T) 


where  the  index  is  inserted  as  a reminder  that  large  AJ  changes  are  pre- 
dicted. In  order  to  extract  the  V-V  rate,  the  experimenter  must  correct  the 
observed  relaxation  rate  for  the  expected  contribution  from  V-R,  T processes. 
In  early  work,  it  was  assumed  that  the  harmonic  oscillator  rule  obeyed;  viz., 

k . = vk.  More  recently,  it  has  been  realized  that  the  v-dependence  is 
v,v-l  1,0  20 

stronger  than  linear;  Kwok  has  suggested  an  exponent  of  2.  3 based  on  ex- 
perimental data  in  several  HF-M  systems.  Wilkins*  numerical  results 
quantitatively  reinforce  Kwok's  suggestion  that  the  corrections  from  V-R,T 
processes  are  large,  and  may,  in  fact,  represent  the  larger  contribution  to 
the  relaxation  rate.  Furthermore,  Wilkins*  finding  that  the  V-V  processes 
involve  changes  of  rotational  state  means  that  one  cannot  use  the  simple, 
rotation-independent  equilibrium  constant  for  the  V-V  process  in  order  to 
convert  the  endothermic  rate  coefficient  into  an  exothermic  one.  For  exam- 
ple, the  equilibrium  constant  for  HF(3)  + HF(0)^  HF(2)  + HF(1),  if  it  is 
assumed  that  species  do  not  change  rotational  states,  is  0.  19.  However, 
Wilkins*  calculations  predict  that  rotational  states  do  change  and,  therefore, 
the  endothermicity  of  the  reaction  is  considerably  smaller  than  if  they  didn't, 
which  results  in  an  equilibrium  constant  of  0.66.  (In  other  words,  the  re- 
action in  the  exothermic  direction  is  only  1.52  times  faster  than  in  the  endo- 
thermic direction,  rather  than  5.2  times  faster.) 


The  calculations  of  Wilkins  have  not,  in  most  cases,  been  directly 
verified  by  clean  experiments.  However,  there  are  considerable  experi- 
mental data  that  strongly  indicate  the  presence  of  high  rotational  energy 
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levels  in  HF  systems,  which  is  what  would  be  expected  if  his  results  are 
correct.  (Wilkins  also  found  that  R-R  relaxation  rates  decrease  noticeably 
with  increasing  J.  Therefore,  although  low  rotational  levels  relax  rapidly, 
high  levels  do  not,  and  this  contributes  further  to  the  perturbations  that  his 
predicted  V-R,  T energy-transfer  processes  would  cause.)  Experiments 
are  being  planned  in  the  Aerophysics  Laboratory  that  will  establish  defini- 
tively the  existence  of  nonequilibrium  populations  of  high  rotational  states 
of  HF,  and,  if  possible,  measure  their  relaxation  rates. 

If  Wilkins'  calculations  are  verified  by  experiments,  there  will  be  two 
important  consequences:  (1)  the  extraction  of  V-V  and  V-R,T  rate  coeffi- 
cients from  experimental  HF  relaxation  data  will  be  even  more  complicated 
than  was  previously  believed,  and  (2)  computer  codes  that  do  not  treat  indi- 
vidual rotational  levels  and  do  not  allow  for  rotational  nonequilibrium  may 
prove  inadequate  for  some  types  of  laser  modeling. 

Kwok  and  Wilkins*  are  working  on  a model  that  incorporates  rotational 
nonequilibrium  but  in  a simplified  manner  so  as  not  to  require  the  full  rota- 
tional manifold  for  every  vibrational  level.  Their  model  approximates  the 
full  rotational  structure  by  incorporating  six  or  seven  discrete  levels  at  high 
J (for  V = 0,  they  retain  J=10  through  16)  and  by  assuming  rotational  equilib- 
rium always  holds  for  the  lower  J levels.  The  justification  for  this  is  that 
the  previous  calculations  of  Wilkins  show  that  rotational  relaxation  is  fast 
for  low  J levels  but  slow  for  upper  J levels. 

Since  the  literature  search  for  Ref.  1 was  completed,  another  major 

experimental  study  on  HF-HF  relaxation  has  been  reported.  Poole  and 
3 1 

Smith  determined  rate  coefficients  for  the  transfer  of  energy  from 
HF(v  = 2.  . , 7)  to  HF(0),  H^i  and  HD  at  room  temperature.  Data  were  ob- 

tained by  monitoring  chemiluminescence  quenching  in  a 10  liter  bulb  into 
which  concentric  flows  of  H^/H/Ar  and  F2/ Ar/ Quencher  were  admitted. 


M.  A.  Kwok  and  R.  L.  Wilkins,  private  communication. 
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Experiments  were  conducted  at  total  pressures  of  3.9  Torr,  with  the  re- 

3 4 

agents  on  the  order  of  1 part  in  10  -10  . First  overtone  emission  was 
monitored  for  v = Z.  . . . 7.  Quenching  of  emission  was  observed  as  a func- 
tion of  pressure  of  added  quenching  gas.  Rate  coefficients  for  deactivation 
by  the  quencher  were  obtained  relative  to  the  sum  of  all  HF{v)  loss  processes 
in  the  absence  of  added  quencher  (i.e.,  due,  principally,  to  spontaneous 
radiation  and,  to  a lesser  extent,  deactivation  by  Ar,  F^,  and  HF 

formed  in  the  reaction).  The  deactivation  rate  coefficients  given  in  Table  5 
were  obtained  for  M=HF,  HD. 


Table  5.  Measured  Deactivation  Rate  Coefficients  for  HF{v) 
by  HF,  H2,  and  D2  at  Room  Temperature 
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^Results  of  Poole  and  Smith  ^3 

Results  of  Kwok  and  Wilkins,  ^ Kwok  and  Cohen,  and  Kwok  (private 
communication) 

^Results  of  Airey  and  Smith^'^ 

Results  of  Osgood,  Sackett,  and  Javan^^ 

^Results  of  Bott^^ 


The  results  up  through  v = 5 obtained  by  Poole  and  Smith  and  by  Kwok  and 
co-workers  are  in  fairly  good  agreement.  As  noted  earlier,  Kwok  has 
maintained  that  HF  V - V relaxation  in  the  presence  of  several  diatomic  and 
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polyatomic  species  obeys  a v dependence  law  of  approximately  v * . The 
results  of  Poole  and  Smith  for  HF  relaxation  in  the  presence  of  HF, 

and  also  fit  a similar  expression  (although  it  is  not  possible  to  choose 

^ 2 2.3 

between,  say,  v and  v * )• 

The  preceding  brief  discussion  should  have  clarified  several  impor- 
tant points : 


1.  There  are  insufficient  experimental  data  with  which  to  evaluate 
individual  rate  coefficients  for  each  V-V  and  V-R,T  process 
involved  in  HF-HF  energy  exchange. 

2.  The  data  that  are  available  are  not  entirely  consistent. 

3.  Theoretical  calculations  provide  a framework  within  which  such 
individual  rate  coefficients  can  be  estimated,  but  the  results  are 
largely  unverified  by  experiment. 

Hence,  we  are  forced  to  combine  experiment  and  theory  in  a somewhat  arbi- 
trary manner  in  order  to  arrive  at  a set  of  rate  coefficients  that  can  be  recom- 
mended for  computer  modeling.  In  Table  6,  we  present  a possible  means  of 
selecting  V-V  and  V-R,  T rate  coefficients  that  is  based  on  the  overall  experi- 
mental deactivation  data  listed  in  Table  5 and  on  the  predictions  of  the  cadcu- 

29 

lations  of  Wilkins.  These  rate  coefficients  are  all  highly  tentative  at  this 
time.  It  should  be  noted  that  they  differ  somewhat  from  the  recommendations 
of  Cohen  and  Bott:^ 

1.  The  V-V  rate  coefficients  of  Ref.  1 were  assumed  to  be  indepen- 
dent of  V.  The  values  given  here  are  inversely  proportional  to  v. 

2.  In  Ref.  1,  it  was  assumed  that  only  single -quantum  V-R,T  de- 
activation occurred,  and  those  rate  coefficients  were  scaled  as 
1:8:28:25:5  for  v=l,2,3,4,  and  5.  Here,  we  permit  multiquantum 
deactivations  to  occur  with  0.5  the  probability  of  single  quantum 
transitions,  and  assume  that  the  total  V-R,  T deactivations  scale 
as  1:6:7:17:30:50. 

The  effect  of  the  change  in  V-V  rates  is  expected  to  be  very  slight  in  pre- 
dicting total  laser  power  output,  although  it  might  affect  noticeably  the  spec- 
tral distribution.  The  effect  of  increased  V-R,  T deactivations  will  mean 
considerably  less  power  predicted  in  all  cases.  Because  the  justification  for 
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Table  6.  Estimation  of  HF(v)  + HF(0)  V-V  and  V-R,  T Rate  Coefficients 
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this  model  is  so  tenuous,  confirmation  by  comparison  with  experimental 
laser  data,  is  recommended. 

For  purposes  of  this  review,  the  best  approach  appears  to  be  to  rely 
on  the  use  of  a model  that  retains  rotational  equilibrium;  there  are  not  suf- 
ficient data  to  justify  the  construction  of  the  full  model  with  separate  rota- 
tional levels.  However,  this  is  done  with  the  recognition  that  continued 
developments  in  laser  science  and  technology  may  render  this  approach 
inadequate  in  another  year  or  two. 

B.  H ATOMS 

Bott  and  Heidner^^  have  extended  their  work  on  HF  relaxation  by 
H and  D atoms  discussed  in  Ref.  1.  They  used  the  same  basic  procedure 
of  laser-induced  fluorescence  measurements  in  a discharge  flow -tube 
apparatus  to  study  rates  of  removal  of  HF(v=l,2,  3)  at  295  K by  H atoms, 
and  the  rate  of  removal  of  HF(3)  by  both  H and  D atoms  at  295,  240,  and 
200  K.  Their  results  for  removal  by  H atoms  at  295  K are  given  in  Table  7 
(from  Ref.  39). 


Table  7. 

Experimental  Removal  Rate  Coefficients 
for  HF(v)  by  H at  T = 295  K 

V 

3 

k(cm  /mole-sec) 

kv/ki 

1 

(1.4  ± 0.4)  xlO^^ 

1.0 

2 

(5  ± 2)  xlO^^ 

3.6 

3 

(6.3  ± 1.5)  xlO^^ 

450 

Removal  rates  for  HF(3)  by  H and  D atoms  at  different  temperatures 
are  given  in  Table  8. 
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Table  8.  Experimental  Removal  Rates  for  HF(3) 
by  H and  D Atoms 


10^^  k,  cm^ /mole -sec 

295  K 

240  K 

200  K 

HF(3)  + H 6.  3 ± 1. 5 

6.  5 ± 0.6 

10. 2 ± 1. 2 

HF(3)  + D 5.9  ± 0.  5 

7.9  ± 1.2 

The  removal  of  HF(v)  by  H atoms  can  be  attributed  to  one  or 

processes,  as  Bott  and  Heidner  have  discussed: 

more  possible 

HF(v)  + H - HF(v')  + H 

v'  <v 

(6-d) 

HF(v)  + H - + F 

(-4) 

The  atom  exchange  process  to  form  is  exothermic  for  v = 3,  but  not  for 
V = 1 or  2 (assuming  formed  in  its  vibrational  ground  state).  The  large 
increase  in  removal  rates  in  going  from  v = 1 to  v = 3 suggests  a contribu- 
tion from  the  atom  exchange  process.  As  Bott  and  Heidner  have  noted,  the 
experimental  rate  coefficient  can  be  interpreted  as  the  sum 


k = ^/  k?,  « /v  + Ak  .. 

exp  V 6(v=3,v  ) -4(v=3) 

where  A is  between  0.7  and  1.0.  From  the  kinetic  and  thermochemical  data 

for  the  reaction  of  F + and  the  known  vibrational  distribution  of  product 

^ 13 

HF,  one  can  calculate  that  k should  be  1.4  x 10  > which  is  about  22% 

of  the  total  removal  rate.  This  suggests  a value  of  between  4.9  and  5.  3 
13 

X 10  for  the  vibrational  deactivation  processes,  i.e.  , for  the  sum 
v H / 

v^^6(v“3  v^)  V = 0,  1,  2.  For  modeling  purposes,  it  is  necessary  to 
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know  the  relative  distribution  among  the  various  values  of  v\  For  example, 
the  effect  on  calculated  laser  performance  of  the  third  level  that  is  being 
deactivated  principally  to  the  second  level  will  be  considerably  different  from 
the  effect  of  deactivation  all  the  way  to  the  zero  level.  Calculations  made  in 
this  Laboratory  indicate  that,  under  some  conditions,  the  assumption  that 
^v  = 3 in  the  relaxation  of  HF(v=3)  by  H atoms  results  in  a prediction  of  ap- 
proximately one-third  of  the  power  predicted  under  the  assumption  of  single- 
quantum deactivation  for  all  HF(v).  In  both  cases,  it  must  be  assumed  that 
the  overall  deactivation  rate  of  HF(v=3)  is  equal  to  the  experimentally 
measured  value  for  HF(3)  deactivation.  Since  the  severe  effects  of  the  multi- 
quantum assumption  seem  to  disagree  with  experience,  we  assume  that  the 
very  large  rate  for  deactivation  of  HF(3)  is  associated  with  single-quantum 
deactivation  only,  assigning  the  rate  by  arbitrarily  selecting  a value  of  0.9 
for  A.  Thus 


H 


k “ k/ / \ "^0.  9k  A*  _o\ 

exp  6(v=3,  V =Z)  -4(v=3) 


H 13  3 

where,  from  figures  given  above,  ^'-2)  “ 5 x 10  cm  /mole-sec. 

Even  after  these  assumptions  have  been  made,  there  is  still  some 
ambiguity  in  the  nature  of  the  deactivation  process.  Deactivation  could  be 
a result  of  a nonreactive,  inelastic  collision,  in  which  the  HF  molecule 
remains  intact  but  loses  vibrational  energy;  or  a result  of  a reactive  colli- 
sion, in  which  the  F atom  is  transferred  to  the  incoming  free  H atom.  Ex- 
periments cannot  distinguish  between  the  two  possibilities,  although  isotopic 
labeling  with  D atoms  should  provide  a good  indication  of  what  is  taking  place 
in  the  H/HF  system.  Wilkins*  trajectory  calculations  can,  in  principle, 
distinguish  between  the  two  processes.  However,  there  is  still  some  uncer- 
tainty in  the  shape  of  the  potential  energy  surface  for  the  H-F-H  system. 
Wilkins  has  studied  the  effect  of  barrier  height  on  the  rate  of  deactivation  of 
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HF(v)  by  H atoms.  He  found,  for  example,  the  approximate  deactivation 
rates  for  HF(v=3)  given  in  Table  9. 


Table  9.  Theoretical  Deactivation  Rates  for  HF(v=3) 
by  H Atoms  for  Different  Assumed  Poten- 
tial Energy  Surface  Barrier  Heights 


Barrier  Height,  kcal/mole 

k,  cm^/mole-sec 

1.  5 

13 

1.5  X 10 

2.  5 

2.  0 X 10^^ 

3.  5 

2.0  X 10^^ 

40.0 

2.  5 X 10^^ 

As  the  barrier  height  increases,  the  contribution  from  reactive  processes 
decreases,  but  the  contribution  from  nonreactive  scattering  increases,  so 
that  the  total  rate  does  not  decrease  monotonically  with  barrier  height. 

No  measurements  have  been  made  on  the  relaxation  of  HF  from  vibra- 
tional levels  higher  than  the  third.  We  assume  that  the  same  large  rate 
coefficient  as  was  found  for  v = 3 -^v  = 2 applies  for  all  v-^v-1  for  v = 3.  We 
assume  further  that  multiquantum  relaxation  occurs  with  the  same  rate  as 
was  found  for  the  deactivation  of  the  v = 2 level.  (The  accuracy  of  this 
estimate  cannot  be  important,  inasmuch  as  the  assumed  rate  coefficient  is 
so  much  smaller  than  that  for  the  single  quantum  process.  ) However,  the 
experimental  disappearance  rate  is  divided  by  2,  which  is  the  number  of 
possible  exit  channels  for  the  case  of  v = 2.  The  assumed  deactivation  coeffi- 
cients are  given  in  Table  10.  The  rates  for  v>3  are  considerably  faster  than 

1 42 

the  earlier  results  of  Kwok  and  Wilkins;  * however,  at  present  we  consider 
these  to  be  the  more  reliable. 
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Table  10.  Suggested  Deactivation  Rate  Coefficients 
for  HF(v)  by  H Atoms 

HF(v)  + H^HF(v')  + H 
12  3 

k = A X 10  exp(-700/RT)cm  /mole-sec 
V v'  A Ref. 


1 0 

2 1,  0 

3 2 


4 

5 

6 

4 

5 

6 


3 

4 

5 


2.  1.0 
3,2,  1,0 
4,  3,2,  1,0 


1.5  Heidner  and  Bott^9 

1.8  Heidner  and  Bott;^*^  with  assump- 

tion that  both  possible  exit  chan- 
nels are  equally  probable. 

360.0  Heidner  and  Bott:^9  experimental 

rate  is  corrected  for  contribution 
from  H_  + F formation  and  then 
assumea  to  represent  single  quan- 
tum deactivation. 


360.0  Assumed  same  as  for  v = 3 


1.8  Assumed  same  as  for  v = 2 


C.  DF 

Recently,  Wilkins*  completed  trajectory  studies  of  the  HF-DF  system, 
in  which  V-V  and  V-R,T  processes  with  either  or  both  species  (HF  and  DF) 
vibrationally  excited  were  taken  into  account.  As  is  the  case  in  many  of  the 
energy-transfer  processes  involving  upper  vibrational  levels  of  HF  in  colli- 
sion with  HF,  these  theoretical  results  have  not  been  verified  by  experimental 
studies  and  should  be  used  with  due  caution.  We  expect,  however,  that  the 


R.  L.  Wilkins,  private  communication. 


-29- 


basic  trends  of  relaxation  probabilities  as  a function  of  varying  vibrational 
level  should  be  reasonably  accurate,  and  we  rely  on  them  in  the  following 
suggested  rate  coefficients. 

The  experimental  data  on  which  the  recommended  value  of  the  V-R,T 
proces  s 


HF(1)  + DF(0)  - HF(0)  + DF(0) 

is  based  were  outlined  in  Ref.  1.  Wilkins  calculated  the  overall  rate 
coefficients  given  in  Table  11  for  HF(v)  + DF(0)  (where  "overall"  means 
summed  over  all  possible  j states). 


Table  11.  Theoretical  Overall  Rate  Coefficients 
for  HF(v)  + DF(0) 


HF(v)  + DF(0)  ^HF(v"  ) + DF(0) 

V 

/ 

V 

k / lO^^cm^/mole-sec  ^/k  /, 

V,  V V V 

10  /mole -sec 

1 

0 

3.6 

3.6 

2 

0 

2.0 

4.  1 

2 

1 

2.  1 

3 

0 

1.5 

3 

1 

1.5 

4.7 

3 

2 

1.7  . 

4 

0 

1.  1 ' 

4 

1 

1.3 

4 

2 

1.  1 

► 

5.0 

4 

3 

1.5 

j 
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A monotonic  increase  with  v,  although  not  large,  is  shown  in  the  last  col- 
umn. For  computational  ease,  and  in  viewof  the  uncertainties  in  the  theo- 
retical results,  we  assume  that  these  values  are  constant  with  increasing  v. 
And  that  the  channels  into  the  various  v*  are  all  equally  probable. 

Column  3 shows  that  the  latter  approximation  is  sufficiently  accurate 
for  most  purposes.  Thus,  the  recommended  rate  coefficient  for  the 
process 

HF(v)  + DF(0)  - HF(v')  + DF(0)  (6-f) 


is 


= - (1.9  X 

O V 


1 2 3 

+ 1.3  X 10  T^) 


The  results  of  Wilkins  for  V-V  energy  transfer  in  the  case  of  HF(v) 
colliding  with  DF(0)  shows  that  the  probability  of  energy  transfer  increases 
significantly  with  v;  his  results  are  given  in  Table  12. 


Table  12. 

Theoretical 
Coefficients 
+ DF(1) 

V-V  Energy  Transfer  Rate 
for  HF(v)  + DF(0)  - HF(v-l) 

V 

k. 

12  3 

10  cm  /mole-sec 

1 

1. 

6 ± 0.  9 

2 

2. 

2 ± 1.0 

3 

3. 

1 ± 1.  2 

4 

4. 

4 ± 1. 4 

Within  the  uncertainty  of  the  calculations,  a linear  dependence  on  v seems 

appropriate.  Hence,  the  recommended  rate  coefficients  are  = 6 xlO^^v 

3 ° 

cm  /mole- sec,  where  reaction  8 is 

HF(v)  ‘+  DF(0)  ^ HF(v-  1)  + DF(  1)  (8) 
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D. 


DEACTIVATION  OF  UPPER  LEVELS  BY 


OTHER  CHAPERONES 
41 

Bott  studied  the  vibrational  relaxation  of  HF(v=l,  2,  3)  in  the  presence 
of  H^f  N2,  and  CO^J  some  results  were  also  reported  for  relaxation  by  O^ 
and  HCl.  The  results  are  summarized  in  Table  13. 


Table  13.  Experimental  Rate  Coefficients  for  Vibrational 
Relaxation  of  HF(v=l,  2,3) 


10 

cm^ /mole  - sec 

V 

^2 

CO2  02 

HCl 

k k /k . 

V V 1 

k 1 

k„/k, 

k 

k /k.  k k /k. 

V 1 V V 1 

k /k, 

V 

V 

V 1 

1 

0.71  1.0 

0.0027 

1.0 

0.  73 

1.0  0.0008  1.0 

0.32  1.0 

2 

1.37  1.9 

0.015 

5.6 

3.7 

5.0 

3 

2.86  4.0 

0.  056 

20.  0 

7.0 

9.6  0.014  16.7 

7.48  23.5 

For  H2>  the  experiments  yield  the  relaxation  rate  in  the  endothermic  di- 
rection, i.e.,  HF(v)  + H2(0)  = HF(v-l)  + H2(l)-  The  figures  given  are  for 
the  reverse,  exothermic  direction,  assuming  no  changes  in  rotational 
energy  levels  of  either  molecule.  See  also  Table  5 for  other  data  for 
relaxation  by  H2* 


V.  ROTATIONAL  RELAXATION 


For  several  years,  spectral  distributions  in  HF  lasers  have  strongly 
suggested  the  presence  of  nonequilibrium  rotational  populations.  The  evi- 
dence was  avoided  as  long  as  possible,  not  so  much  because  it  appeared 
physically  unreasonable  as  because  it  appeared  computationally  intractable 
on  an  affordable  timescale. 

Results  of  the  first  direct  measurements  of  rotational  relaxation  in  HF 

indicated  that  the  process  was  extremely  rapid,  with  rate  coefficients  on  the 
14  3 43-45 

order  of  10  cm  /mole-sec.  The  extraction  of  these  rate  coefficients 

from  experimental  data  required  some  assumptions  about  the  various  pro- 
cesses that  were  taking  place  in  the  medium.  In  general,  processes  ex- 
changing rotational  energy  with  translational,  or  transferring  rotational 
energy  from  one  molecule  to  another,  were  considered;  processes  that  ex- 
change rotational  energy  with  vibrational  energy  were  neglected.  An  example 
of  the  latter  is 


HF(v=l,  J=3)  + HF(v=0,  J = 16)  + M 

Recent  calculations  by  Wilkins  have  suggested  that  such  processes  can  be 
very  rapid,  which  makes  doubtful  the  rate  coefficients  that  had  been  deduced 
from  the  experimental  data.  Furthermore,  both  experimental  and  theoretical 
results  indicate  that  rotational  relaxation  rates  decrease  with  increasing  J. 
Therefore,  even  if  rotational  energy  transfer  were  to  occur  with  nearly  gas - 
kinetic  frequency  for  low  J,  it  could  be  orders  of  magnitude  slower  for  higher 

* 3}C 

J levels  of  interest.  Calculations  by  Kerber  and  by  Hough  indicate  that, 
even  if  rotational  energy  transfer  were  to  occur  on  every  gas-kinetic  colli- 
sion, that  would  still  not  be  fast  enough  to  maintain  rotational  equilibrium  in 
a lasing  mediiim. 


R.  L.  Wilkins,  private  communication. 
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The  picture,  then,  appears  rather  unclear  now,  although  it  does  appear 
clear  that  rotational  nonequilibrium  is  an  issue  that  must  be  taken  into  con- 
sideration. However,  it  is  not  obvious  that  a model  needs  to  be  constructed 
that  treats  each  rotational  level  separately.  Kwok  and  Wilkins  have  been 
investigating  a kinetic  model  that  sets  Idwer  rotational  levels  in  equilibrium 
(on  the  assumption  that  energy  transfer  among  them  is  very  rapid),  but 
treats  upper  levels  discretely.  The  degree  of  success  of  this  approach  has 
not  yet  been  fully  determined. 

As  noted  earlier,  the  recommended  rate  coefficients  given  in  the 
Appendix  are  based  on  data  interpreted  with  the  assumption  of  rotational 
equilibrium.  The  extent  to  which  this  assumption  detracts  from  the  useful- 
ness of  the  computer  modeling  remains  to  be  seen. 
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APPENDIX 

RECOMMENDED  RATE  COEFFICIENTS  FOR  - F^ 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
Y pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 

of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
I sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 

i cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  ele ctro- optic s ; communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

' Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 

matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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